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Engineers  (OCE) ,  U.  S.  Army,  16  July  1980,  at  the  request  of  the  U.  S. 
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Simmons,  Chief  of  the  Hydraulics  Laboratory,  J.  L.  Grace,  Jr.,  Chief  of 
the  Hydraulic  Structures  Division,  and  Dr.  D.  R.  Smith,  Chief  of  the 
Reservoir  Water  Quality  Branch  (Physical).  The  study  was  conducted  by 
Mr.  J.  P.  Holland  with  some  assistance  from  Mr.  B.  Loftis.  The  report, 
prepared  by  Mr.  Holland  with  assistance  from  Ms.  B.  R.  Turner,  was  re¬ 
viewed  by  Mr.  Grace  and  Dr.  Smith. 

Commanders  and  Directors  of  WES  during  the  conduct  of  this  study 
and  the  preparation  and  publication  of  this  report  were  COL  Nelson  P. 
Conover,  CE,  and  COL  Tilford  C.  Creel,  CE.  Technical  Director  was 
Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units 

verted  to  metric  (SI) 

Mult  ip ly 

of  measurement  used 

units  as  follows: 

_ By _ 

in  this  report  can  be  con- 

To  Obtain 

acre-feet 

1233.482 

cubic  metres 

acres 

4046.856 

square  metres 

Btu  (International  Table) 

1055.056 

joules 

cubic  feet  per  second 

0.02831685 

cubic  metres  per  second 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins* 

feet 

0.3048 

metres 

feet  per  second 

0.3048 

metres  per  second 

feet  per  second  per 
second 

0.3048 

metre:.  oer  second  per  second 

miles  (U.  S.  statute) 

1.609344 

kilometres 

pounds  (mass)  per 
cubic  foot 

16.01846 

kilograms  per  cubic  metre 

square  feet 

0.09290304 

square  metres 

square  miles  (U.  S. 
statute) 

2.589988 

square  kilometres 

*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  read¬ 
ings,  use  the  following  formula:  C  =  (5/9) (F  -  32).  To  obtain  Kelvin 
(K)  readings,  use:  K  =  (5/9)(F  -  32)  +  273.15. 
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EFFECTS  OF  STORAGE  REALLOCATION  ON  THERMAL 


CHARACTERISTICS  OF  COWANF.SQUE  LAKE,  PENNSYLVANIA 
Numerical  Model  Investigation 

PART  I:  INTRODUCTION 

Background  and  Purpose 

1.  Growing  consumptive  use  of  water  in  the  Susquehanna  River 
Basin  has  resulted  in  several  users  expressing  an  interest  in  purchasing 
water  supply  storage  in  Cowanesque  Lake,  a  project  which  is  nearing  the 
filling  stage  of  construction.  The  lake,  as  originally  planned,  has  no 
water  supply  storage.  The  U.  S.  Army  Engineer  District,  Baltimore  (NAB), 
is  examining  the  feasibility  of  reallocating  some  flood-control  storage 
of  Cowanesque  Lake  to  water  supply.  This  reallocation  would  result  in  a 
substantial  increase  in  normal  pool  (30  to  40  ft*)  elevation  which  may 
alter  the  thermal  structure  of  the  lake.  Consequently,  present  selec¬ 
tive  withdrawal  capabilities  may  prove  inadequate.  The  purpose,  then, 

of  this  study  is  twofold.  The  first  is  the  investigation  of  the  temper¬ 
ature  structure  expected  in  Cowanesque  Lake  after  the  normal  pool  is  in¬ 
creased  following  reallocation  of  some  flood-control  storage  to  water 
supply  storage  within  the  existing  project.  The  second  purpose  of  the 
study  is  the  location  of  additional  selective  withdrawal  intakes  that 
would  allow  operation  of  the  system  to  satisfy  downstream  warmwater  tem¬ 
perature  objectives. 

Project  Description 

2.  The  Cowanesque  Lake  project  is  located  on  the  Cowanesque  River 
approximately  2.2  miles  upstream  of  its  confluence  with  the  Tioga  River 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (SI)  units  is  presented  on  page  4. 
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at  Lawrencevil  le,  Pennsylvania,  and  approxima^t  I.-'  12  miles  south  of 
Corning,  New  York,  as  shown  in  Figure  1.  Cowanesque  Dam,  which  is  of 
rolled  earth  and  rock-fill  construction,  is  3100  ft  in  length  with  a 
crown  151  ft  above  the  streambed.  At  the  dam  there  is  a  drainage  area 
of  293  square  miles.  The  project  was  designed  to  impound  water  primarily 
for  flood  control.  Recreational  use  was  also  authorized. 

3.  When  the  lake  is  at  the  elevation  of  the  spillway  crest,  el 
1117.0,*  it  will  extend  approximately  8.0  miles  upstream  and  will  have  a 
surface  area  of  approximately  2060  acres.  For  the  initially  designed 
normal  recreation  pool  (el  1045.0),  the  lake  has  a  maximum  depth  of  45  ft 
and  extends  4.2  miles  upstream  with  a  surface  area  of  410  acres. 

4.  Releases  from  the  lake  will  be  made  through  an  outlet  works 
with  an  outlet  tunnel  whose  cross-sectional  area  is  139  sq  ft  and  over  a 
400-ft-long  emergency  spillway.  The  outlet  works  tunnel  is  excavated 
through  rock  undek.  the  right  abutment  of  the  dam.  The  outlet  works 
presently  consist  of  an  intake  control  tower  with  provisions  for  making 
selective  withdrawal  and  flood  releases.  Flood  releases  through  the  out¬ 
let  works  will  be  controlled  by  two  6-f t-wide  by  14-ft-high,  hydrauli¬ 
cally  operated,  slide-type  service  gates.  The  selective  withdrawal  re¬ 
leases  will  pass  through  two  identical  wet  wells,  each  housing  two  5.0- 
by  5.0-ft  intakes,  one  at  center-line  el  1014.5  and  one  at  center-line 

el  1037.0.  All  four  intakes  discharge  into  the  front  face  of  their  re¬ 
spective  wet  well.  The  flow  from  each  wet  well  will  be  controlled  by 
2-ft-wide  by  4-ft-high,  hydraulically  operated  slide  gates  located  in 
separate  passages  on  both  sides  of  the  flood-control  passages,  down¬ 
stream  of  the  service  gates.  These  slide  gates  are  at  center-line 
el  1014.0.  Each  intake  will  pass  a  maximum  of  250  cfs  at  normal  pool 
(el  1045.0). 

Approach 

5.  This  study  was  accomplished  with  the  use  of  numerical  modeling. 
The  approach  involved  the  selection  of  ^veral  study  years  and  simulation 

*  All  elevations  (el)  cited  herein  are  in  feet  referred  to  the  National 
Geodetic  Vertical  Datum  (NGVD). 
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(><•  ?!  ion  «ind  virinirv  maps 


of  lake  operation  for  each  of  these  years.  Study  years  selected  had  com 
binations  of  strearaflow  quantities  and  air  temperatures  that  could 
create  varying  conditions  of  thermal  stratification.  The  input  data  re¬ 
quired  by  the  model  were  lake  inflows  and  outflows,  inflow  stream  temper 
atures,  meteorological  conditions  for  each  of  the  study  years,  geometry 
of  the  lake,  and  geometry  of  the  intake  structure. 

G.  The  heat  transfer  into  and  out  of  the  lake  was  evaluated  and 
the  heat  was  distributed  within  the  lake.  A  heat  budget  was  computed 
continuously  throughout  each  simulation  period.  An  objective  tempera¬ 
ture  was  specified  for  each  day,  and  an  operating  scheme  was  determined 
to  best  achieve  that  objective.  The  operation  scheme  for  any  day  was 
the  combination  of  open  ports  that  minimized  the  difference  between  the 
objective  downstream  temperature  and  the  predicted  release  temperature. 
The  output  from  the  computations  included  comparisons  between  objective 
and  release  temperatures  in  graphical  form  throughout  the  simulation 
period  as  well  as  tabular  summaries  for  each  day  and  plotted  profiles  of 
temperature  within  the  lake  at  specified  times  of  the  year. 
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PART  II:  DISCUSSION  OF  NUMERICAL  METHODOLOGY 


7.  Location  of  additional  selective  withdrawal  intakes  for  the 
Cowanesque  Lake  thermal  study  was  accomplished  through  use  of  a  hybrid 
numerical  approach  which  involved  the  coupling  of  a  thermal  simulation 
model  with  a  mathematical  optimizer.  An  overview  of  the  components  of 
this  numerical  approach  appears  below. 

Thermal  Model  Description 

8.  The  downstream  release  temperatures  and  the  in-lake  tempera¬ 
ture  characteristics  for  Cowanesque  Lake  were  predicted  using  a  thermal 
simulation  model.  The  model  (WESTEX)*  used  in  conjunction  with  this  in¬ 
vestigation  was  developed  at  the  U.  S.  Army  Engineer  Waterways  Experi¬ 
ment  Station  (WES)  based  on  the  results  of  Clay  and  Fruh  (1970),  Edinger 
and  Geyer  (1965),  Dake  and  Harleraan  (1966),  and  Bohan  and  Grace  (1973). 

9.  The  WESTEX  model  provides  a  procedure  for  examining  the  bal¬ 
ance  of  thermal  energy  imposed  on  an  impoundment.  This  energy  balance 
and  lake  hydrodynamic  phenomena  are  used  to  map  vertical  profiles  of 
temperature  in  the  time  domain.  The  model  includes  computational  methods 
for  simulating  heat  transfer  at  the  air-water  interface,  heat  advection 
due  to  inflow  and  outflow,  ami  internal  dispersion  of  thermal  energy. 
These  computational  components  require  various  hydrological  and  meteoro¬ 
logical  data  that  are  explained  in  the  following  paragraphs.  A  thorough 
discussion  of  the  WESTEX  model  appears  in  Appendix  A. 

Development  of  Thermal  Model  Inputs 

10.  As  stated  in  the  previous  paragraph,  the  WESTEX  thermal  model 
requires  input  data  for  lake  inflows  and  outflows,  inflow  stream  tempera¬ 
ture,  and  various  meteorological  conditions.  Further,  calibration  of 
the  model  requires  a  determination  of  the  appropriate  surface  exchange 


*  Bruce  Loftis,  '*WESTEX-A  Reservoir  Heat  Budget  Model"  (first  draft), 
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and  internal  mixing  coefficients.  This  calibration  procedure  and  the 
various  inputs  to  the  thermal  model  are  described  below  for  the 
Cowanesque  study. 


Study  Years 

11.  The  years  studied  in  this  investigation  were  specified  by  NAB 
as  follows: 

£.  1958  -  high-flow  year. 

b_.  1964  -  low-flow  year. 

£.  1966  -  low-flow  year. 

d_.  1967  -  average-flow  year. 

Simulations  were  run  from  March  through  October  for  each  of  these  years. 
In  general,  it  is  this  period  in  which  density  stratification  occurs. 

Meteorology 


12.  Meteorological  data  from  Williamsport,  Pennsylvania,  Class  A 
Weather  Station  were  obtained  from  the  National  Oceanic  and  Atmospheric 
Administration  for  use  in  this  study.  Williamsport  is  approximately 

40  miles  south  of  the  Cowanesque  project.  The  required  data  consist  of 
dry  bulb  and  dewpoint  temperatures,  wind  speed,  and  cloud  cover.  Daily 
average  values  were  computed  for  each  and  used  to  determine  equilibrium 
temperatures,  surface  heat  exchange  coefficients,  and  daily  average  solar 
radiation  quantities  for  the  period  of  record.  NAB  supplied  these  three 
parameters  for  1958,  1966,  and  1967;  1964  values  were  computed  at  WES. 

Hydrology 

13.  Hydrologic  routings  were  conducted  by  NAB  to  determine  ex¬ 
pected  daily  inflow  and  outflow  quantities.  These  are  shown  in 

Plates  1A  and  IB  for  each  of  the  study  years.  Outflows  for  the  low-flow 
years  of  1964  and  1966  included  augmented  low-flow  releases  in  order  to 
meet  water  supply  constraints  downstream;  1958  and  1967  required  no 
augmented  releases  for  water  supply. 
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Lake  Geometry 


14.  The  area-capacity  curve  for  Cowanesque  Lake  is  shown  in  Fig¬ 
ure  2.  This  curve  and  other  data  describing  the  location  and  design  of 
the  intake  structure  were  furnished  by  NAB. 

AREA,  HUNDREDS  OF  ACRES 


0  20  40  60  80  100  120  140  160  180  200  220 

CAPACITY,  THOUSANDS  OF  ACRE-FT 

Figure  2.  Area-capacity  curve 

Natural  Downstream  Temperature 

15.  In  order  to  establish  a  natural  downstream  temperature  objec¬ 
tive,  NAB  performed  a  least  squares  curve  fit  to  the  observed  stream 
temperatures  for  all  years  of  record.  The  resulting  equation  was 

0  =  A  sin  (Bt  +  C)  +  D  (1) 
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where 


0  =  natural  stream  temperature  objective  on  day  t,  °F 
t  =  calendar  day  of  year 
A  =  -21.685 
B  =  0.0172 
C  =  1.0908 
D  =  52.579 


Model  Calibration 


16.  The  WESTEX  model  requires  the  determination  of  coefficients 
of  surface  heat  exchange  distribution  and  internal  mixing.  For  Cowan- 
esque  Lake  these  coefficients  were  determined  by  conducting  simulations 
with  Cowanesque  hydrologic  and  meteorologic  data.  Coefficients  were  ad¬ 
justed  and  simulation  was  repeated  until  the  predicted  temperature  pro¬ 
files  corresponded  in  shape  and  range  to  those  observed  by  NAB  personnel 
in  North  Central  Pennsylvania  lakes  of  similar  size,  depth,  and  flow 
magnitude.  The  following  coefficients  were  determined  from  the  analysis: 


P  =  0.5 
A.  =  0.17 
ctj  =  0.30 
<*2  =  0.20 


where 


P  =  percentage  of  incoming  shortwave  radiation  absorbed  in  sur¬ 
face  layer 

A  =  light  extinction  coefficient,  ft  ^ 

Ofj  =  mixing  coefficient  at  surface 
a2  =  mixing  coefficient  at  bottom 


The  A  ,  p  ,  and  values  were  found  to  agree  very  well  with  the 

values  used  in  the  thermal  simulation  model  study  of  Kinzua  Lake  (Dortch 
1981),  a  project  in  the  same  geographical  area. 
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Formulation  of  P< art  Location  Procedure 

17.  The  principal  objective  of  the  Cowanesque  Lake  thermal  study 
was  to  determine  the  location  of  additional  selective  withdrawal  intakes 
in  order  to  meet  a  specified  downstream  warmwater  objective.  To  achieve 
this  objective,  an  optimization  algorithm  was  coupled  with  the  WESTEX 
model  to  predict  the  optimum  port  elevations  for  a  given  set  of  con¬ 
straints.  These  "optimum”  elevations  are  the  set  of  port  locations 
whose  operation  yielded  the  minimum  difference  between  predicted  release 
and  objective  temperatures  over  the  simulated  period.  The  optimum  com¬ 
bination  of  port  elevations  was  defined  as  that  which  would  result  in 
the  minimum  objective  function.  The  objective  function  selected  basi¬ 
cally  consisted  of  the  sum  of  squared  deviations  between  the  computed 
release  and  objective  temperatures.  The  components  of  the  port  location 
procedure  are  described  in  the  following  sections. 

Port  Location  Procedure 

18.  For  each  individual  operating  condition,  the  location  of  addi¬ 
tional  selective  withdrawal  intakes  was  accomplished  through  a  sequential 
process.  This  process  involved  estimation  of  the  optimum  intake  config¬ 
uration  for  downstream  temperature  maintenance,  prediction  of  downstream 
release  characteristics  resulting  from  operation  of  this  intake  config¬ 
uration,  and  computation  of  an  objective  function  (described  in  para¬ 
graph  19).  This  process  was  repeated  until  a  minimum  value  for  the  ob¬ 
jective  function  was  obtained,  at  which  point  the  estimated  location  of 
the  selective  withdrawal  intakes  was  defined  as  the  "optimum"  location 
for  maintenance  of  downstream  temperature  objectives.  Although  this 
procedure  could  be  accomplished  manually,  an  optimization  routine  was 
used  to  locate  an  estimate  to  the  "optLmum"  selective  withdrawal  config¬ 
uration  for  each  operating  condition.  The  port  location  procedure  used 
in  this  study  is  described  in  detail  in  Appendix  B. 

Development  of  Objective  Function 

19.  The  principal  purpose  of  this  study  was  to  determine  the 
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location  of  additional  selective  withdrawal  intakes  to  meet  specified 
downstream  warmwater  temperature  objectives.  An  objective  function  is  a 
scalar  index  that  provides  an  indication  of  how  good  one  possible  deci¬ 
sion  (i.e.,  location  estimate)  is  in  meeting  these  given  objectives. 
Minimization  of  this  objective  function  yields  the  optimal  location  of 
these  additional  ports  for  release  temperature  control. 

20.  The  objective  function  used  for  Cowanesque  Lake  is  the  sum  of 
several  terms.  Two  terms  are  weighted  portions  of  the  sum  of  squares  of 
deviations  between  predicted  release  and  target  temperatures.  The  prin¬ 
cipal  term  represents  the  contribution  of  the  squared  temperature  devia¬ 
tion  on  days  when  the  absolute  value  of  the  deviation  (predicted  release 
temperature  from  target  temperature)  is  greater  than  2.78°C  (5°F).  This 
portion  has  physical  significance  in  that  the  State  of  Pennsylvania  re¬ 
quites  releases  from  Cowanesque  Lake  to  be  plus  or  minus  2.78°C  (5°F) 
from  the  natural  stream  temperature  (which  is  the  study  objective).  A 
number  of  other  objective  function  terms,  which  were  added  to  enhance 
numerical  uniqueness  and  reflect  physical  constraints,  are  described  in 
Appendix  B. 
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PART  III:  NUMERICAL  SIMULATIONS 


21.  A  broad  range  of  conditions  was  simulated  during  the 
Cowanesque  study.  Simulations  were  conducted  for  normal  pool  elevations 
1075.0,  1080.0,  and  1085.0  for  each  of  four  study  years  (1958,  1964, 

1966,  and  1967).  The  three  pools  represent  probable  normal  pools  follow¬ 
ing  reallocation  of  a  portion  of  flood-control  storage  to  water  supply 

s  torage . 

22.  In  addition  to  these  simulated  years  and  pools,  NAB  presently 
has  three  possible  schemes  for  the  operation  of  the  Cowanesque  project 
for  water  supply  releases.  These  three  schemes,  or  "scenarios,"  repre¬ 
sent  allocation  of  water  supply  releases  for  differing  consumers  down¬ 
stream.  The  three  scenarios  are  characterized  by  the  following  flow 
targets  downstream: 

Scenario  1-50  cfs  target  at  Wilkes-Barre,  Pa. 

Scenario  2-45  cfs  target  at  Harrisburg,  Pa. 

Scenario  3-50  cfs  target  at  Wilkes-Barre,  Pa. 

45  cfs  target  at  Harrisburg,  Pa. 

These  three  scenarios  were  simulated  for  each  pool  for  the  1964  study 
year.  Only  Scenario  2  was  simulated  for  1966,  as  requested  by  NAB. 

Study  years  1958  and  1967  (average-  and  high-flow  years,  respectively) 
required  no  water  supply  releases. 

23.  Simulations  of  each  combination  of  pool,  year,  and  scenario 
were  made  for  each  of  the  three  normal  pools  with  operation  of  the  ex¬ 
ist.  ng  selective  withdrawal  intake  configuration.  Since  these  simula¬ 
tions  showed  temperature  objectives  would  not  be  met,  additional  with¬ 
drawal  ports  were  studied.  Analysis  by  NAB  personnel  indicated  that  up 
to  three  additional  ports  could  be  placed  in  each  of  two  wet  wells  at 
the  Cowanesque  project.  Thus,  up  to  six  separate  port  locations  could 
be  sited  in  addition  to  those  existing  in  order  to  meet  specified  con¬ 
straints.  Analogous  simulations  incorporating  additions  to  the  existing 
selective  withdrawal  system  were  then  run. 

24.  The  geometry  of  these  additional  ports  was  another  decision 
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variable.  Ports  sized  the  same  as  those  presently  in  place  at 
Cowanesque,  5  ft  wide  by  5  ft  high  (hereafter  referred  to  as  small), 
were  used  for  most  of  this  study.  Ports  6  ft  wide  by  7  ft  high  (here¬ 
after  referred  to  as  large),  which  allowed  additional  selective  with¬ 
drawal  capacity  beyond  that  of  the  small  ports,  were  also  analyzed. 

25.  Discussion  of  these  various  simulations  appears  in  PART  IV. 


PART  IV:  DISCUSSION  OF  RESULTS  FOR  SIMULATION-SPECIFIC 

CONDITIONS 

26.  As  discussed  in  PART  III,  a  broad  range  of  conditions  was 
simulated  during  the  Cowanesque  study.  Combinations  of  differing  pool 
elevations,  study  years,  port  geometry  and  number,  and  release  scenarios 
were  simulated.  For  each  of  these  combinations,  an  "optimum"  intake  con¬ 
figuration  for  that  specific  combination  was  found.  This  PART  discusses 
results  from  simulation  of  these  "specific"  port  configurations.  Of 
primary  importance  in  this  discussion  are  the  number  and  location  of  ad¬ 
ditional  intakes  required  for  each  simulation  to  meet  the  specified  down¬ 
stream  temperature  objective.  These  results  may  be  divided  into  the 
following  areas: 

a.  No  additional  ports  added  to  the  existing  system. 

b. .  One  additional  level  of  ports  added  to  the  existing 

system. 

£.  Two  additional  levels  of  ports  added  to  the  existing 
system. 

Coupled  with  the  discussion  of  the  number  and  location  of  additional 
ports  is  the  presentation  of  simulation  results  for  specific  conditions 
relating  the  effects  of  both  release  "scenario"  strategies  and  larger 
port  geometry  on  downstream  temperature  maintenance.  Selective  with¬ 
drawal  system  flexibility  requirements  for  the  specific  conditions  simu¬ 
lated  are  discussed  at  the  end  of  this  PART. 

No  Additional  Ports 

27.  Operation  of  the  present  selective  withdrawal  configuration 
(hereafter  referred  to  as  0-levels)  for  each  pool,  year,  and  release 
scenario  combination  was  simulated.  Results  of  these  simulations,  shown 
graphically  as  plots  of  daily  release  temperature  in  Plates  2A-2C,  show 
the  inadequacy  of  the  present  system  to  meet  the  specified  downstream  tem¬ 
perature  objective.  The  two  solid,  smooth  curves  on  each  plot  trace  the 
boundary  of  the  ±2.78°C  (±5°F)  band  around  the  natural  stream  temperature 
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objective.  The  release  temperatures  resulting  from  operation  of  the 
present  selective  withdrawal  system  were  generally  outside  this  ±2.78°C 
band  from  late  spring  to  late  summer  for  each  year.  The  number  of  days 
the  daily  release  temperature  violated  this  band  appears  in  Table  1  for 
each  simulation.  A  minimum  of  100  violation  days  was  computed  for  the 
present  system. 


One  Additional  Level 


28.  The  inadequacy  of  the  present  selective  withdrawal  system  to 
meet  downstream  temperature  objectives  prompted  investigation  of  addi¬ 
tional  port  elevations  in  orde1*  to  meet  the  warmwater  temperature  objec¬ 
tive.  Initial  attempts  to  locate  additional  ports  showed  that  the 
siting  of  one  port  level  (two  ports  at  the  same  elevation,  one  port  per 
wet  well)  was  efficient  in  meeting  downstream  objectives.  Results  from 
simulation  of  the  Cowanesque  selective  withdrawal  system  with  one  addi¬ 
tional  level  of  5-ft-wide  by  5-ft-high  ports  (referred  to  as  1-level 
small  ports  in  this  section)  are  presented  graphically  in  Plates  2A-2C. 

A  summary  of  the  center-line*  elevations  of  the  additional  level  intakes, 
as  well  as  the  number  of  days  temperatures  released  through  this  system 
violated  the  ±2.78°C  objective  band,  appears  in  Table  2.  As  indicated 
in  Plates  2A-2C,  the  addition  of  one  level  of  small  ports  placed  at  the 
"optimum  locations"  of  Table  2  provided  vast  improvement  in  the  ability 
of  the  selective  withdrawal  system  to  meet  the  specified  downstream  tem¬ 
perature  objective  when  compared  with  operation  cf  the  existing  system. 
Release  temperatures  are  rarely  outside  the  objective  temperature  band 
following  addition  of  a  single  level  of  small  ports.  When  the  predicted 
temperature  is  outside  the  objective  band,  it  is  usually  cooler  than  the 
objective.  As  shown  in  Tab1^  2,  the  average  difference  between  release 
and  objective  temperatures  is  consistently  negative.  The  duration  of 
these  occurrences  is  generally  small,  typically  on  the  order  of  10  days 


*  All  port  elevations  in  this  and  subsequent  sections  are  port  center- 
line  elevations. 


18 


or  less.  The  magnitude  and  duration  of  these  release  temperature  devia¬ 
tions  are  presented  graphically  in  Plates  3A-3C.  Plotted  on  these 
pi-  es  are  the  algebraic  differences  between  the  predicted  release  tem¬ 
perature  and  the  objective  temperature  for  each  simulated  day. 

29.  The  discussion  presented  above  was  for  the  results  of  simula¬ 
tion  of  one  additional  level  of  ports  located  at  the  elevations  given  in 
Table  2.  These  ’'optimum"  port  elevations  should  not,  however,  be  con¬ 
strued  as  "absolute"  with  respect  to  maintenance  of  the  given  downstream 
temperature  objective.  For  the  set  of  conditions  simulated,  a  range  of 
center-line  elevations  exists  over  which  little  change  in  predicted  re¬ 
lease  temperature  occurs.  This  is  indicated  in  plots  of  the  objective 
function  (described  in  paragraph  19)  as  a  function  of  port  level  eleva¬ 
tion  as  shown  in  Plates  4A-4C.  For  each  condition  a  range  of  center- 
line  intake  elevations  exists  over  which  little  change  in  the  objective 
function  value  around  the  minimum  value  is  noted.  In  general,  for  years 
with  no  water  supply  releases  and  little  pool  fluctuation  (1958  and 
1967),  this  flat  range  in  the  function  is  found  for  port  elevations  from 
near  the  normal  pool  water  surface  to  approximately  10  ft  below  the 
normal  pool  (approximately  the  minimum  thickness  of  the  epilimnion,  the 
upper,  warmwater  layer,  for  the  simulation  period).  For  1964  and  1966, 
years  with  water  supply  releases,  this  range  runs  from  approximately 

10  ft  below  to  approximately  20  ft  below  the  elevation  of  the  normal 
pool  water-surface  elevation.  Thus,  for  each  simulation  there  is  ap¬ 
proximately  a  minimum  10-ft  range  in  which  the  center-line  elevation  of 
the  additional  selective  withdrawal  intake  level  may  be  placed  without  a 
loss  in  system  ability  to  meet  downstream  temperature  objectives. 

Two  Additional  Levels 

30.  For  the  simulated  conditions,  a  single  level  of  ports  strate¬ 
gically  located  in  the  pool  resulted  in  a  selective  withdrawal  design 
that  would  usually  release  water  within  the  desired  temperature  band. 
When  the  release  temperature  was  outside  the  objective  band  it  was  con¬ 
sistently  too  cold.  In  general,  to  improve  temperature  maintenance  with 
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a  second  level  of  additional  ports,  it  would  be  necessary  to  site  the 
second  level  above  the  elevation  of  the  optimum  single  level  locations 
(shown  in  Table  2)  so  that  more  warm  water  could  be  withdrawn.  This  ap¬ 
proach  was  employed  in  the  evaluation  of  the  addition  of  a  second  level 
port  to  meet  the  desired  warmwater  objective. 

31.  For  all  pool  elevations  investigated  for  1958  ar,d  1967  a>  1 
for  the  el  1075.0  pool  of  1964  (all  three  scenarios),  no  improvement  in 
meeting  the  temperature  objective  was  indicated  by  the  addition  of  two 
port  level  elevations  as  compared  with  one.  An  analysis  of  the  meteoro¬ 
logical  and  hydrological  data  for  these  simulations  helps  to  explain 
this  conclusion.  In  the  spring,  when  both  warm  temperatures  and  warm 
intermediate  inflows  would  normally  promote  the  blending  of  water  from 
the  upper  layers,  the  high  outflows  simulated  significantly  exceeded  the 
capacity  of  the  selective  withdrawal  system,  thereby  requiring  the  re- 

ase  of  cold  bottom  waters  from  the  floodgates.  The  location  of  a 
greater  number  of  additional  intakes  was  precluded  by  the  effects  of 
thermal  blockage  which  reduces  or  negates  flow  control  of  individual  in¬ 
takes  which  are  operating  simultaneously  in  the  same  wet  well.  Further, 
during  the  summer  the  simulated  outflows  did  not  exceed  the  selective 
withdrawal  capacity  of  a  single  level  of  intakes.  As  a  result,  only  one 
additional  level  of  intakes  was  required  to  produce  release  temperatures 
consistent  with  downstream  summer  temperature  objectives  for  these  simu¬ 
lations. 

32.  A  lack  of  improvement  in  temperature  maintenance  with  the  ad¬ 
dition  of  a  second  port  level  a?  described  above  was  also  noted  for  the 
1075.0  pool  elevation  for  1966.  The  lack  of  improvement,  however,  was 
not  related  to  those  causes  described  above.  The  specific  meteorological 
and  hydrological  conditions  of  1966,  coupled  with,  the  1075.0  pool  eleva¬ 
tion,  resulted  in  little  difference  between  the  temperature  near  the 
water  surface  and  optimum  single  port  level.  The  low  spring  outflows 
were  generally  less  than  the  selective  withdrawal  capacity  of  a  single 
level  of  intake.  Thus,  no  improvement  in  release  quality  was  predicted 
through  the  operation  of  a  second  additional  port. 
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33.  For  the  simulations  of  i 9 6 A  (all  three  scenarios)  and  1966  in 
which  the  pool  elevations  were  1080.0  and  1085.0,  there  was  an  improve¬ 
ment  in  the  ability  of  the  selective  withdrawal  system  to  meet  the  down¬ 
stream  temperature  objective  through  the  use  of  two  additional  levels  of 
intakes  as  compared  with  a  single  level.  Improvement  occurred  during 
the  midspring  period  (mid-May  to  early  June).  During  this  period,  be¬ 
fore  the  onset  of  water  supply  releases,  the  pool  remains  steady  at 

el  1085.0  or  el  1080.0  for  the  respective  study  year.  In  each  case, 
meteorological  and  hydrological  conditions  produced  a  thermal  profile 
which  was  generally  1  to  2.5°C  wanner  near  the  surface  than  at  the  eleva¬ 
tion  of  the  optimum  single  level  location.  For  respective  simulated 
conditions,  the  optimum  single  port  location  was  approximately  10  to 
20  ft  below  the  normal  pool  elevation  to  compensate  for  summer  pool 
drawdown.  For  the  flow  conditions  simulated,  ports  located  10  to  20  ft 
(center-line  elevation)  below  the  normal  pool  surface  did  not  withdraw 
adequate  quantities  of  warm  epilimnetic  waters  during  the  spring  period, 
while  ports  located  near  the  normal  pool  surface  could  withdraw  these 
waters. 

Scenario  Differences 

34.  NAB  designated  three  water  supply  release  strategies,  or 
scenarios,  as  described  in  paragraph  22.  Of  the  years  studied,  all 
three  scenarios  were  used  only  in  1964.  Plots  of  both  daily  release 
temperature  and  deviation  of  this  temperature  from  the  objective  for 
1964  appear  in  Plates  2A-2C  and  3A-3C,  respectively.  Examination  of 
these  plates  revealed  no  significant  difference  in  the  thermal  character¬ 
istics  of  Cowanesque  Lake  releases  as  a  result  of  the  differing  1964  re¬ 
lease  strategies  of  Scenarios  1,  2,  and  3.  Further,  no  difference  in 
the  location  of  additional  upper-pool  port  levels,  or  in  the  objective 
function  values,  was  noted  among  the  scenarios.  Because  differences  be¬ 
tween  the  scenarios  were  considered  minor,  only  Scenario  3  was  run  for 
1964  during  the  remainder  of  the  sti  'v. 
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Large  Port  Geometry 


35.  Examination  cf  the  release  temperature  deviations  incurred 
through  the  use  of  the  optimum  port  configuration  for  a  specific  simu¬ 
lation  condition  revealed  a  strong  correlation  between  large  release 
temperature  deviations  and  the  release  of  flows  in  excess  of  the  selec¬ 
tive  withdrawal  system  capacity.  To  minimize  these  deviations,  the 
capabilities  of  the  Cowanesque  selective  withdrawal  system  should  be 
maximized  within  the  structural  constraints  of  the  Cowanesque  intake 
structure.  As  an  example,  when  the  pool  is  at  el  1085.0,  the  wet  wells 
can  pass  a  maximum  combined  flow  of  700  cfs.  However,  due  to  frictional 
and  entrance  losses,  a  5-f t-wide  by  5-ft-high  port  of  the  type  presently 
installed  at  Cowanesque  will  pass  only  250  cfs,  resulting  in  500  cfs 
maximum  for  combined  flow  through  the  selective  withdrawal  system  with 
these  ports.  Thus,  only  approximately  71  percent  of  the  potential  re¬ 
lease  capacity  of  the  Cowanesque  selective  withdrawal  system  would  be 
utilized  with  5-f t-wide  by  5-ft-high  ports  for  pool  elevation  of  1085.0. 

36.  To  maximize  the  selective  withdrawal  capabilities  of  the  sys¬ 
tem,  simulations  were  made  with  a  single  level  of  additional  ports  6  ft 
wide  by  7  f t  high.  Examination  of  the  Cowanesque  outlet  works  and  dis¬ 
cussion  with  NAB  personnel  Indicated  that  this  port  geometry  was  a 
aropcr  size  to  facilitate  both  the  desired  flow  capacity  and  the  hydrau¬ 
lic  and  structural  constraints  of  the  Cowanesque  wet  wells.  The  ports 
were  placed  at  the  optimum  center-line  location  found  for  the  5-ft-wide 
by  5-ft-high  single  level  as  given  by  Table  2.  These  configurations 
were  then  simulated  for  each  study  year  and  pool.  The  selective  with¬ 
drawal  maximum  outflow  was  set  at  700,  670,  and  640  cfs  for  the 

el  1085.0,  1080.0,  and  1075.0  pools,  respectively.  Daily  release  tem¬ 
perature  plots  for  these  large  ports  are  shown  in  Plates  5A-5C.  Objec¬ 
tive  function  values  and  the  number  of  days  of  temperature  violations 
for  large  port  operation  are  given  in  Table  3. 

37.  Differences  in  the  in-lake  thermal  distribution  as  a  result 
of  the  use  of  large  rather  than  small  ports  were  also  noted.  Plates  6A- 
6X  show  the  in  situ  temperature  profiles  expected  at  Cowanesque  Lake  for 
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representative  simulations  of  both  large  and  small  ports.  These  profiles 
show  that  during  1958,  1964,  and  1967,  a  1  to  2°C  cooling  of  the  hypolim- 
nion  (bottom,  cold  water  layer)  was  predicted  through  the  use  of  large 
ports  when  compared  with  similar  small  port  operation.  This  results  from 
the  release  of  larger  quantities  of  epilimnion  waters  in  late  spring  by 
the  6-ft-wide  by  7-ft-high  ports  than  by  the  smaller  ports.  Larger  epi¬ 
limnion  release,  in  turn,  resulted  in  the  release  of  smaller  amounts  of 
hypolimnion  waters  through  the  floodgate.  Thus,  cold  hypolimnion  waters 
that  would  have  been  released  and  subsequently  replaced  by  wanner  waters 
from  higher  in  the  pool  remain  within  the  pool  when  the  large  ports  are 
used.  This  cooling  trend  remains  in  effect  throughout  the  majority  of 
the  simulation  period,  resulting  in  warmer  bottom  temperatures  through 
use  of  the  5-ft-wide  by  5-ft-high  ports  than  through  use  of  the  larger 
counterpart.  Surface  temperatures,  however,  are  consistently  the  same 
for  both  port  types.  The  1966  study  year  had  no  large  spring  flows,  as 
observed  in  the  other  study  years,  which  necessitated  hypolimnion  re¬ 
leases.  Thus,  the  operational  schemes  of  both  the  6-ft-wide  by  7-ft- 
high  ports  and  the  5-ft-wide  by  5-ft-high  ports  were  identical  for  the 
period,  and  the  cooling  effect  described  abjve  was  not  observed  for  1966. 

38.  The  impact  to  aquatic  biota  as  a  result  of  the  in-lake  cool¬ 
ing  experienced  by  large  port  use  relative  to  small  port  use  is  beyond 
the  scope  of  this  work.  If  this  in-lake  cooling  can  be  accepted,  the  use 
of  larger  port  geometries  of  a  size  capable  of  meeting  the  maximum  selec¬ 
tive  withdrawal  capacity  and  the  physical  constraints  of  the  Cowanesque 
wet  wells  seems  warranted.  The  improvements  in  meeting  the  downstream 
obj  ctive  through  use  of  the  large  ports  may  outweigh  predicted  degrada¬ 
tions  in  the  lake.  Further,  the  larger  vertical  dimension  of  the  6-ft- 
wide  by  7-ft-high  ports  would  allow  greater  flexibility  in  both  fully  and 
partially  submerged  operation.  This  flexibility  would  be  of  great  impor¬ 
tance  during  periods  with  large  pool  fluctuations  such  as  those  of  1964. 

Multiport  Flexibility 

39.  Although  one  level  of  additional  selective  withdrawal  intakes 
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is  generally  adequate  to  meet  downstream  temperature  objectives  for  each 
of  the  specific  combinations  considered,  there  are  conditions  when  a 
single  level  located  at  the  optimum  elevation  does  not  effectively  main¬ 
tain  temperature  objectives.  In  instances  where  the  optimum  single  level 
of  ports  was  located  10  to  20  ft  below  the  normal  pool  elevation,  such  as 
1964  or  1966,  test  results  showed  that  an  additional  port  level  located 
high  in  the  epilimnion  would  result  in  improvement  to  spring  release  tem¬ 
perature  maintenance  beyond  that  of  a  single  level  lower  in  the  pool. 
Further,  due  to  the  relatively  small  incremental  volumes  of  Cowanesque 
Lake,  large  net  inflows  will  result  in  a  substantial  increase  in  pool 
elevation.  Routings  for  1964  indicate  that  a  net  inflow  of  13,100  cfs 
(on  day  65)  would  result  in  an  increase  of  approximately  20  ft  in  the 
pool.  The  location  of  a  level  of  ports  near  the  normal  pool  water  sur¬ 
face  in  addition  to  those  located  at  a  lower  optimum  single  elevation 
would  facilitate  the  meeting  of  downstream  temperature  at  increased  pool 
elevations. 

40.  Instances  may  also  exist  wherein  an  optimum  single  level  of 
ports  located  near  the  normal  pool  elevation  (as  in  1958  and  1967)  would 
not  effectively  maintain  release  temperature  objectives.  Large  pool 
drawdown,  such  as  that  for  1964,  would  render  a  level  of  ports  located 
near  the  normal  pool  water  surface  inoperative  over  a  large  portion  of 
the  simulation  period.  This  action  would  result  in  very  inefficient  re¬ 
lease  of  epilimnion  waters  due  to  the  depth  of  the  highest  existing  port 
level  (el  1037.0),  thereby  increasing  the  magnitude  and  frequency  of  tem¬ 
perature  violations.  Further,  flexibility  would  also  be  needed  during 
periods  of  filling  when  a  single  level  of  ports  high  in  the  pool  would 
preclude  epilimnetic  release  by  selective  withdrawal  prior  to  port  sub¬ 
mergence.  Intermediate  ports  to  this  upper  level  of  intakes  and  those 
existing  are  requisite  to  system  flexibility  during  filling. 
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PART  V:  OPTIMUM  MULTIPORT  CONFIGURATION 


41.  In  PART  IV,  results  were  presented  which  showed  that  a  single 
additional  set  of  intakes  located  strategically  in  the  pool  would  gener¬ 
ally  provide  adequate  downstream  temperature  maintenance  for  the  specific 
conditions  simulated.  However,  results  also  indicated  that  conditions 
can  exist  which  necessitate  the  location  of  two  sets  of  additional  in¬ 
takes  rather  than  a  single  additional  set.  If  two  sets  of  intakes  are 
added  to  the  present  system,  selective  withdrawal  system  flexibility  for 
a  eiv?’i  year  and  normal  pool  elevation  could  then  be  substantially  en¬ 
hanced  through  use  of  multiple  port  elevations  with  the  following  port 
scheme : 

£.  Two  large  ports  (one  level)  located  high  in  the  epilimnion 
in  a  given  specific  range  for  the  normal  pool  and  year 
simulated. 

_b.  One  large  port  located  at  the  elevation  of  the  optimum 

1964  single  intake  level  for  the  given  normal  pool  (shown 
in  Table  2). 

c.  One  large  port  located  10  ft  below  the  optimum  1964  eleva¬ 
tion  of  b_. 

d^.  Small  ports  of  the  existing  system. 

The  lower  additional  ports  (described  in  _b  and  £  above)  would  be  used 
only  during  periods  of  extreme  pool  drawdown.  In  general,  the  releases 
simulated  during  these  drawdown  periods  were  less  than  the  maximum  selec¬ 
tive  withdrawal  flow  (350  cfs)  through  a  single  large  port.  Therefore, 
single  large  ports  located  at  each  of  these  two  elevations  for  a  speci¬ 
fied  pool  would  be  adequate  to  pass  the  simulated  outflows. 

42.  The  placement  of  the  additional  ports  in  the  manner  described 
above  would  result  in  the  following: 

a.  System  flexibility  would  be  increased  for  all  normal  pools 
simulated  for  1964.  Placement  of  the  highest  additional 
intakes  in  the  epilimnion  of  the  pool  would  help  reduce 
spring  release  temperature  degradations.  Placement  of  an 
additional  intake  near  the  optimum  single  elevation  for 
the  given  normal  pool  would  result,  by  definition,  in  the 
minimum  release  temperature  deviation  from  the  objective 
for  the  bulk  of  the  simulation  period. 
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t>.  Placement  of  an  additional  intake  near  the  optimum  1964 
single  elevation  would  allow  the  highest  additional  level 
of  intakes  for  the  1080.0  and  1078.0  pools  of  1966  to  be 
placed  higher  in  the  pool  than  their  optimum  single  level 
elevations.  Results  indicated  that  the  highest  addi¬ 
tional  level  could  be  placed  2  and  4.5  ft  above  the  opti¬ 
mum  1966  single  level  elevation  (center-line  el  1071.4 
and  el  1066.0  for  pool  elevations  1080.0  and  1075.0,  re¬ 
spectively)  with  no  loss  in  temperature  maintenance  for 
these  pool  elevations,  respectively. 

£.  Epilimnetic  releases  of  at  least  350  cfs  could  be  achieved 
for  all  pool  elevations  between  the  normal  operating  pool 
elevation  and  an  elevation  25  ft  below.  Thus,  a  25-ft 
drawdown  (which  was  exceeded  only  once  in  this  studv) 
could  be  incurred  with  little  loss  in  selective  withdrawal 
efficiency  for  outflows  less  than  350  cfs. 


Pool-Specific  Optimum  Multiport  Configuration 


43.  Locating  four  additional  ports  in  the  manner  described  in 
paragraph  41  results  in  changing  the  range  of  elevations  over  which  the 
highest  level  of  additional  ports  could  be  located.  This  range  is  given 
for  both  the  addition  of  a  single  level  of  ports  and  the  addition  of 
four  ports  (as  u«_ov.ribed  in  paragraph  41)  in  Table  4.  As  stated  in  para¬ 
graph  42,  the  range  is  extended  for  the  1030.0  and  1075.0  pools  of  1966 
and  all  pools  for  1964.  A  consequence  of  this  broadened  range  for  the 
two  low-flow  years  is  that  port  locations  can  be  selected  which  are  com¬ 
mon  to  all  conditions  simulated  for  a  given  normal  pool  elevation. 

Table  4  shows  that  the  optimum  four-port  configuration  which  is  common 
to  all  study  years  for  each  of  the  three  normal  pools  consists  of  the 
fol  lowing: 


b. 


1085.0  pool. 

Two  intakes  located  between  el  1075.0  and  el  1077.0. 

One  port  located  at  el  1067.0. 

One  port  located  at  el  1057.0. 

Intakes  of  existing  system  (two  ports  at  el  1037.0,  two 
at  el  1014.5). 

1080.0  pool. 

Two  intakes  located  between  el  1070.0  and  el  1073.5. 
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One  port  located  at  el  1063.0. 

One  port  located  at  el  1053.0. 

Intakes  of  existing  system. 
c_.  1075.0  pool. 

Two  intakes  located  between  el  1065.0  and  el  1071.5. 

One  port  located  at  el  1062.0. 

One  port  located  at  el  1052.0. 

Intakes  of  existing  system. 

If  the  upper  level  of  intakes  is  located  at  el  1065.0  for  the  1075.0 
pool,  the  port  at  el  1062.0  could  be  placed  at  el  1052.0  or  removed  with 
little  loss  of  selective  withdrawal  efficiency  for  conditions  simulated. 
It  should  be  noted  that  all  ports  except  those  existing  were  assumed  to 
be  6  ft  wide  by  7  ft  high. 

Optimum  Port  Configuration  for  All  Pool  Elevations 

44.  If  the  normal  operating  pool  elevation  is  known  to  be  a  spe¬ 
cific  pool  from  among  the  three  studied,  the  selective  withdrawal  con¬ 
figuration  outlined  in  the  previous  paragraph  for  the  particular  pool 
would  be  an  optimum  scheme  for  temperature  maintenance  downstream.  If 
this  normal  pool  elevation  is  not  well  defined,  or  if  the  normal  pool 
elevation  is  later  changed,  a  selective  withdrawal  system  designed  for  a 
specific  pool  elevation  could  prove  inadequate.  An  example  of  such  ac¬ 
tion  would  be  the  change  in  normal  pool  elevation  from  el  1085.0  to 
el  L075.0  which  would  leave  the  highest  level  of  ports  inoperative  over 
a  large  portion  of  the  simulation  period.  Thus,  the  possibility  of  one 
single  configuration  of  intakes  which  could  be  operated  optimally  for 
temperature  maintenance  for  all  of  the  conditions  studied  was  investi¬ 
gated.  Determination  of  this  configuration  (which  was  to  consist  of 
four  large  ports  in  addition  to  those  existing)  was  begun  by  selecting 
an  optimum  elevation  for  the  highest  additional  level  of  ports  which  was 
common  to  all  pool/year  conditions.  This  elevation  was  determined  by 
examining  the  range  of  optimum  center-line  elevations  given  in  Table  4 
for  the  highest  additional  ports  for  each  pool/year  condition  with  only 
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a  single  level  of  ports  added.  Results  for  the  addition  of  a  single 
level  of  ports  rather  than  those  for  the  addition  of  four  ports  were  con¬ 
sidered  because  the  latter  results  were  based  on  locations  for  the  third 
and  fourth  additional  ports  which  were  pool-specific  and  thus  contrary 
to  this  effort.  Examination  of  Table  4  shows  that  a  common  range  of  op¬ 
timum  center-line  elvations  for  the  highest  additional  ports  extends 
from  el  1070.0  to  el  1071.4  for  all  conditions  simulated  except  the 
1958/1085.0  pool  combination  and  the  1075.0  pool  combinations  of  1964  and 
1966.  The  upper  end  of  the  range,  el  1071.4,  was  chosen  for  the  eleva¬ 
tion  of  the  highest  additional  intakes.  Two  additional  intakes,  one  each 
at  el  1062.0  and  el  1052.0,  were  then  sited  for  additional  flexibility. 
This  configuration  is  also  one  of  the  optimum  four-port  configurations 
for  the  1075.0  pool  (described  in  paragraph  43)  and,  as  such,  will  be 
quite  adequate  at  meeting  downstream  objectives  for  the  1075.0  pool 
combinations  of  1964  and  1966.  Representative  results,  Table  5  and 
Plates  7A-7C  and  8A-8X,  show  that  with  the  exception  of  the  1958/1085.0 
pool,  this  system  is  capable  of  meeting  downstream  temperature  objec¬ 
tives  at  least  as  well  as  the  optimum  port  configuration  for  any  of  the 
specific  year/pool  elevation  combinations.  For  hydrological  and  meteoro¬ 
logical  conditions  similar  to  those  simulated  in  this  study,  this  single 
configuration  of  ports  is  generally  adequate  to  meet  downstream  tempera¬ 
ture  constraints. 

45.  The  1958/1085.0  pool  results  indicate  the  limits  of  tl  is  con¬ 
figuration.  Location  of  the  upper  port  level  at  el  1071.4  did  not  allow 
adequate  withdrawal  of  epilimnetic  waters  for  1958  (the  high  inflow, 
high  outflow  year).  This  occurred  because  the  majority  of  the  epilim- 
nion  rose  above  the  effective  withdrawal  limit  of  the  1071.4  port  eleva¬ 
tion.  However,  it  should  be  noted  that  no  selective  withdrawal  config¬ 
uration  simulated  was  highly  efficient  in  meeting  downstream  temperature 
objectives  for  the  high  inflow/outflow  conditions  of  1958.  The  upper 
intake  level  at  el  1071.4  is  within  the  epilimnion  for  all  other  simu¬ 
lated  conditions,  and  no  temperature  deviation  in  excess  of  that  expected 
with  a  "pool-specific"  configuration  should  be  expected.  Thus,  with  the 
exception  of  high-flow  years,  a  single  combination  of  port  locations  may 


be  used  to  meet  the  desired  temperature  objective  for  all  conditions 
similiar  to  those  modeled. 

Selective  Withdrawal  Capacity 

46.  The  capacity  of  the  selective  withdrawal  system  at  Cowanesque 
Lake  is  700  cfs  at  el  1085.0  normal  pool.  While  this  capacity  is  gener¬ 
ally  adequate,  hydrologic  events  during  the  spring  of  three  of  the  study 
years  (1966  excluded)  resulted  in  an  extended  period  wherein  the  capacity 
of  the  selective  withdrawal  system  was  inadequate  to  meet  the  downstream 
objei  ive.  NAB  has  indicated  that  it  plans  to  maintain  a  maximum  normal 
pool  at  el  1085.0  and  that  storage  above  this  elevation  will  be  quickly 
released.  If  the  releases  from  the  top  ports  are  at  a  maximum,  700  cfs, 
any  additional  flow  must  be  released  from  the  floodgate.  This  flow, 
which  is  extremely  cold  in  the  spring,  may  cause  a  temperature  deviation 
in  the  spring  of  up  to  2°C  over  a  one-day  time-step  for  a  floodgate  flow 
as  little  as  300  cfs.  Any  upgrading  of  the  selective  withdrawal  capacity 
would  be  beneficial.  The  Cowanesque  system  seems  especially  sensitive 
to  bottom  releases  in  the  spring,  and  this  increase  in  selective  with¬ 
drawal  capacity  may  merit  consideration. 
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PART  VI:  CONCLUSIONS  AND  RECOMMENDATIONS 


47.  A  numerical  simulation  model  in  connection  with  a  mathemati¬ 
cal  optimizer  has  been  used  to  locate  additional  selective  withd’awal 
intakes  at  Cowanesque  Lake.  These  additional  intakes  will  be  needed  in 
order  to  meet  a  specified  downstream  warmwater  temperature  objective  fol¬ 
lowing  increases  in  normal  pool  elevation  as  a  result  of  proposed  reallo¬ 
cation  of  some  flood-control  storage  to  water  supply  storage.  Operation 
of  the  present  selective  withdrawal  system  showed  the  system  to  be  in¬ 
adequate  in  meeting  temperature  objectives  under  increased  pool  eleva¬ 
tion  conditions.  Simulations  indicated  that  a  single  level  of  addi¬ 
tional  ports,  whose  center-line  elevations  were  located  from  near  the 
normal  pool  water  surface  to  approximately  10  ft  below  that  water  sur¬ 
face,  was  adequate  (except  during  extremely  high  Inflow  periods)  to  meet 
downstream  objectives  for  study  years  1967  and  1958  (average-  and  high- 
flow  years,  respectively).  The  upper  bound  of  this  10-ft  region  was  ap¬ 
proximately  10  ft  below  the  normal  pool  water  surface  for  study  years 
1964  and  1966  (low-flow  years). 

48.  Simulation  indicated  that  the  majority  of  temperature  devia¬ 
tions  observed  occurred  when  the  outflow  was  in  excess  of  the  selective 
withdrawal  capacity.  The  use  of  larger  ports  (6  ft  wide  by  7  ft  high) 
allowed  an  increase  in  the  selective  withdrawal  capacity  from  Its  pres¬ 
ent  500  cfs  (as  constrained  by  the  5-ft-wide  by  5-ft-high  existing  ports) 
to  700  cfs.  Either  improved  or  unchanged  temperature  deviations  were 
generally  noted  through  use  of  the  larger  port  size  for  the  additional 
intakes  when  compared  with  the  5-ft-wide  by  5-ft-high  ports. 

49.  Results  from  various  simulations  indicated  a  need  for  a  sec¬ 
ond  set  of  additional  intakes.  These  Intakes  would  allow  a  vastly  more 
flexible  selective  withdrawal  system  than  a  single  level  of  additional 
intakes  would  afford.  This  flexibility  would  be  required  during  periods 
of  highly  fluctuating  pool  elevations.  Further  flexibility  would  be 
needed  during  the  filling  process,  when  a  single  additional  port  level 
located  high  above  the  bottom  would  not  provide  the  necessary  selective 
withdrawal  c  ipabilities  for  pool  elevations  between  its  elevation  and 
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the  elevation  of  the  highest  existing  ports  (el  1037.0). 

30.  It  is  recommended  that  two  additional  sets  of  6-ft-wide  by 
7-ft-high  selective  withdrawal  intakes  be  added  to  the  existing  ports  in 
the  Cowanesque  Dam  outlet  works.  If  the  normal  operating  pool  elevation 
is  known  from  the  three  elevations  studied,  it  is  recommended  that  two 
upper  level  intakes  be  placed  with  their  center  lines  approximately  3  to 
10  ft  below  the  elevation  of  that  water  surface.  Placement  of  additional 
6-ft-wide  by  7-ft-high  ports  at  the  optimum  single  level  location  corre¬ 
sponding  to  1964  conditions  with  the  normal  pool  under  consideration 
would  allow  a  common  range  of  optimum  upper  level  elevations  to  exist 
for  all  simulations  of  a  given  pool.  If  the  hydrologic  conditions  simu¬ 
lated  in  this  study  are  representative  (spring  storm  events,  small  sum¬ 
mer  outflows  during  periods  of  lake  drawdown),  a  single  port  at  the  1964 
optimum  location  would  suffice.  The  fourth  additional  6-ft-wide  by  7-ft- 
high  port  could  then  be  placed  with  a  center  line  10  ft  below  the  1064.0 
single-level  optimum  elevation  to  extend  system  flexibility. 

51.  If  the  normal  operating  pool  cannot  be  defined  at  this  time, 
or  if  the  normal  pool  might  be  changed  later,  a  different  intake  config¬ 
uration  is  suggested.  This  configuration  is: 

<i.  Two  large  ports  (one  level)  at  el  1071.4 
b^.  Single  large  port  at  el  1062.0 
£.  Single  large  port  at  el  1052.0 
d.  Ports  of  the  existing  system 

where  all  port  elevations  are  center-line  elevations.  The  upper  intake 
level  may  not  withdraw  sufficient  warm  waters  from  the  epilimnion  during 
high-flow  years.  High-flow  years  in  combination  with  a  1085.0  pool,  or 
other  hydrological  and  meteorological  factors  not  specifically  simulated 
in  this  study,  could  negate  the  effectiveness  of  this  alternate  config¬ 
uration.  If  the  simulated  hydrological  and  meteorological  conditions 
are  representative,  then  this  system  is  quite  adequate  for  all  pools  and 
years  simulated  except  the  combination  of  high  inflow  years  with  a  normal 
pool  at  el  1085.0. 

52.  Regardless  of  the  intake  configuration,  it  is  recommended  that 
all  additional  intakes  located  at  Cowanesque  Lake  be  of  sufficient 
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dimension  to  pass  the  maximum  selective  withdrawal  capacity  attainable 
for  the  pools  studied,  700  cfs.  This  selective  withdrawal  capacity  may 
still  be  inadequate  for  spring  outflows  of  as  little  as  300  cfs  in  ex¬ 
cess  of  the  present  selective  withdrawal  capacity. 

53.  The  simulations  for  all  pools  of  the  1958  study  year  (high 
flow)  showed  a  marked  inability  to  meet  downstream  temperature  objec¬ 
tives  in  the  spring  season.  A  period  of  spring  temperature  violations 
was  also  noted  for  1967  (average  flow).  These  violations  are  a  result 
of  a  combination  of  factors.  The  selective  withdrawal  capacity  at 
Cowanesque  Lake,  as  discussed  in  the  previous  paragraph,  is  often  too 
small  to  release  the  required  flow.  Secondly,  there  were  periods  in  mid 
spring,  such  as  in  1958,  when  no  large  quantities  of  warm  water  existed 
in  the  Cowanesque  pool.  Both  of  these  factors  are  contributors  to  re¬ 
lease  temperature  violations  during  periods  of  high  spring  outflow.  Out 
flow  not  withdrawn  through  upper  pool  outlets  is  released  through  the 
floodgates.  Release  of  this  very  cold  hypolimnetic  water  drives  the  re¬ 
lease  temperature  well  below  that  of  the  upper  pool  temperatures.  The 
1967  violations  were  of  this  nature.  When  extremely  high  inflow  and  out 
flow  conditions  are  coupled  with  extremely  cold  inflow  temperatures,  the 
ability  of  the  present  Cowanesque  selective  withdrawal  system  to  meet  a 
specified  release  temperature  is  further  lessened.  The  severity  and 
duration  of  the  release  temperature  deviations  observed  for  1958  are  an 
example. 
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Table  1 

Results  of  Simulation  of  Existing 
Selective  Withdrawal  System 


Average 


Normal 

Pool 

Elevation 
Year  ft  NGVD 

Scenario 

Objective 

Function 

Value 

Number  of 
Days  of 
Band 

Violation 

Maximum 

Band 

Violation 

°C 

Dif  ference 
from 

Objective 

°C 

1967 

1085 

NA* 

7,140 

141 

9.80 

-4.39 

1958 

1085 

NA 

w  » - -  “ 

148 

9.16 

-4.18 

1966 

1085 

2 

12,700 

161 

11.92 

-5.74 

1964 

1085 

1 

6,610 

139 

9.52 

-3.93 

1964 

1085 

2 

7,129 

144 

9.52 

-4.19 

1964 

1085 

3 

6,550 

136 

9.52 

-3.83 

1967 

1080 

NA 

4,500 

135 

8.03 

-3.53 

1958 

1080 

NA 

7,139 

142 

10.83 

-4.34 

1966 

1080 

2 

9,727 

151 

10.59 

-4.95 

1964 

1080 

1 

4,737 

119 

8.36 

-3.24 

1964 

1080 

o 

5,120 

126 

8.36 

-3.46 

1964 

1080 

3 

4,729 

118 

8.36 

-3.16 

1967 

1075 

NA 

3,184 

115 

7.09 

-2.88 

1958 

1075 

NA 

5,554 

147 

10.20 

-3.93 

1966 

1075 

2 

7,051 

142 

9.50 

-4.15 

1964 

1075 

1 

3,077 

100 

7.15 

-2.53 

1964 

1075 

2 

3,076 

100 

7.15 

-2.53 

1964 

1075 

3 

3,077 

100 

7.15 

-2.48 

*  NA  denotes  not  applicable.  No  water  supply  releases  for  these  years 


*  Port  configuration  consisted  of  two  ports  at  the  elevation  of  the  single  optimum,  plus  ports  of  the 
existing  system  (two  ports  at  el  1037.0,  two  ports  at  el  1014.5). 

**  NA  denotes  not  applicable.  No  water  supply  release  for  these  years. 
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Table  5 


anH  All  Ports  of  rh»  Existing  System  _ 


Year 


1967 

1958 

1966 
1964 

1967 
1958 

1966 
1964 

1967 
1958 
1966 
1964 


Normal 

Pool 

Elevation 
ft  NGVD 

1085 

1085 

1085 

1085 

1080 

1080 

1080 

1080 

1075 

1075 

1075 

1075 


Scenario 

NAI¬ 

NA 

2 

3 

NA 

NA 

2 

3 

NA 

NA 

2 

3 


Objective 

Function 

Value 

135 

712 

163 

16 

103 

459 

105 

19 

80 

447 

89 

12 


Number  of 
Days  of 
Band 

Violation 

4 

46 

14 

1 

4 

27 
10 

1 

-4 

28 
8 
1 


Maximum 

Band 

Violation 

°C 


6.72 

7.88 

4.13 

3.90 

5.82 

7.60 

3.80 

4.29 


4.77 

7.58 

3.83 

3.43 


Average 

Difference 

from 

Objective 

°C 


-0.55 

-1.51 

-0.88 

-0.21 

-0.36 

-1.20 

-0.59 

-0.19 

-0.36 

-1.07 

-0.57 

-0.18 


ft. 


*  All  port  elevations  are  center-line^^vationSo  ^  ^ 

“  applicable .  -ter  supply  releases  for  these  years 
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Plate  2B.  Daily  release  temperatures  for  the  existing  system 
(0-levels)  and  existing  system  plus  one  additional  level  of 
small  ports  (1-level),  pool  el  1080.0 
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Plate  3A.  Deviation  of  predicted  daily  release  temperature 
from  objective  temperature  with  one  additional  level  of 
small  ports,  pool  el  1085.0 
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SCENARIO  2 
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ICIGHT  OF  AOOITIOMAL  POST  LEVEL  ABOVE  BOTTOM.  FT 


ICIGHT  OF  AOOITIOMAL  PORT  LEVEL  ABOVE  BOTTOM.  FT 


Plate  4A.  Objective  function  value  as  a  function  of  port  height 
above  lake  bottom  (el  1000.0)  with  one  additional  level  of  small 

ports,  pool  el  1085.0 
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SCENARIO  2 
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tCOMT  OF  ADDITIONAL  PORT  LEVEL  ABOVE  BOTTOM.  FT 


(OVNT  OF  ADDITIONAL  PORT  LEVEL  ABOVE  BOTTOM,  FT 
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SCENARIO  3 


tCIGHT  OF  ADDITIONAL  PONT  LEVEL  A30VE  0OTTW.  FT 


PlAt-e  AB  Objective  function  value  as  a  function  of  port  height 

t*.  hot  tom  (el  1000.0)  with  one  aOditlonnl  level  of  -.11 
ports,  pool  el  lOoU.U 
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Plate  IB.  Release  temperatures  for  one  additional  level  of  ports 
by  5-  by  5-ft  (small)  ports  and  one  additional  level  of  ports  of 
6-  by  7-ft  (large)  ports,  pool  el  1080.0 
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Plate  6B.  Temperature  profiles,  pool  el  1085.0,  Jul-Oct  1958,  with  one 
additional  level  of  large  and  small  ports 


HCIQHT  mOM  BOTTOM,  FKT 


TEMrewrnjRE,°c 


TEHPDWTURE,  “C 


Plate  6F.  Temperature  profiles,  pool  el  1085.0,  Jul-Oct  1966,  with  one 
additional  level  of  large  and  small  ports 
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Plate  6P.  Temperature  profiles,  pool  el  1080.0,  Jul-Oct  1967,  with  one 
additional  level  of  large  and  small  ports 
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Plate  6Q.  Temperature  profiles,  pool  el  1075.0,  Mar-Jun  1958,  with  one 
additional  level  of  large  and  small  ports 
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Plate  6R.  Temperature  profiles,  pool  el  1075.0,  Jul-Oct  1958,  with  one 
additional  level  of  large  and  small  ports 
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Plate  6T.  Temperature  profiles,  pool  el  1075 .0,  Jnl  uct 

additional  level  of  large  and  small  ports 
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Plate  6W.  Temperature  profiles,  pool  el  1075.0,  Mar-Jun  1967,  with  one 
additional  level  of  large  and  small  ports 
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Plate  8B.  Temperature  profiles  with  single  optimum  port 
configuration,  pool  el  1085.0,  Jul-Oct  1958 
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Plate  8C.  Temperature  profiles  with  single  optimum  port 
configuration,  pool  el  1085.0,  Mar-Jun  1964 
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Plate  8E.  Temperature  profiles  with  single  optimum  port 
configuration,  pool  el  1085.0,  Mar-Jun  1966 
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Plate  8F.  Temperature  profiles  with  single  optimum  port 
configuration,  pool  el  1085.0,  Jul-Oct  1966 
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Plate  8T.  Temperature  profiles  with  single  optimum  port 
configuration,  pool  el  1080.0,  Mar-Jun  1958 
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Plate  81,.  Temperature  profiles  with  single  optimum  port 
configuration,  pool  el  1080.0,  Jul-Oct  1964 
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Plate 


8N.  Temperature  profiles  with  single  optimum  port 
configuration,  pool  el  1080.0,  Jul-Oct  1966 
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8R.  Temperature  profiles  with 
configuration,  pool  el  1075.0, 


single  optimum  port 
Jul-Oct  1958 
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APPENDIX  A:  DISCUSSION  OF  WESTEX  THERMAL  MODEL 


Simulation  Model  Description 

1.  The  downstream  release  characteristics  and  internal  structure 
of  temperature  within  a  reservoir  are  predicted  with  a  numerical  simula¬ 
tion  model.  The  model,  hereinafter  identified  as  WESTEX,  was  developed 
by  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  based  on 
results  of  Clay  and  Fruh  (1970),  Edinger  and  Geyer  (1965),  Dake  and 
Harleman  (1966),  and  Bohan  and  Grace  (1973).* 

Introduction 


2.  The  reservoir  is  conceptualized  as  a  number  of  homogeneous 
horizontal  layers  stacked  vertically,  and  the  neat  sources  and  sinks  to 
a  general  layer  are  represented  as  shown  in  Figure  Al.  The  solution  for 
the  temperature  history  of  a  general  layer  is  obtained  by  solving  a  con¬ 
servation  of  mass  and  energy  equation.  The  governing  equation  is: 


90.  0.Q-  0  Q  .  a  /  90t\  .  9(Q  0.)  ,  au 

L  _  o^o  13  L\  1  ^v  LI  9H 

at  aaz  aaz  a  az  ^  az  j  *  a  '  az  pcpA  az 

where 

0^  =  temperature  of  layer,  °F 
t  =  time,  days 
0.  =  inflow  temperature,  °F 

^  ''i 

Q.  =  flow  rate  into  layer,  ft  /day 

1  ,  2 
A  =  horizontal  cross-sectional  area,  ft 

AZ  =  layer  thickness,  ft 

0  =  outflow  temperature,  °F 

°  3 

Qq  =  outflow  rate,  ft  /day 

Z  =  elevation,  ft 

2 

k  =  vertical  diffusion  coefficient,  ft  /day 


(Al) 


*  See  REFERENCES  at  end  of  main  text. 


Al 


CONTROL  VOLUME 
REPRESENTATION 


Figure  Al.  Typical  layer  in  one-dimensional  description 


3 

Q  =  net  vertical  flow  into  or  out  of  layer,  ft  /day 

v  3 

p  =  density  of  water,  lb/ft 

=  specific  heat  of  water,  Btu/lb/°F 
H  =  external  heat  source,  Btu/ft/day 
Appropriate  boundary  conditions  roust  be  supplied  for  inflow  and  outflow 
rates,  inflow  temperatures,  and  surface  heat  exchange  at  the  air-water 
interface.  Solution  of  Equation  A1  for  each  layer  through  time  yields 
the  dynamic  vertical  temperature  distribution  of  the  reservoir.  Funda¬ 
mental  assumptions  and  various  processes  addressed  to  solve  these  types 
of  problems  will  be  discussed  in  the  following  sections. 

Fundamental  Assumptions 

3.  Reservoir  hydrodynamic  phenomena  and  .a  thermal  energy  balance 
are  used  to  predict  temperature  profiles  and  release  temperatures  in  the 
time  domain.  The  model  includes  computational  methods  for  simulating 
heat  transfer  at  tne  air-water  interface,  advective  heat  due  to  inflow 
and  outflow,  and  the  internal  diffusion  of  thermal  energy.  The  model  is 
conceptually  one-dimensional  based  on  the  division  of  the  impoundment 
into  discrete  horizontal  layers  of  uniform  thickness.  Assumptions  in¬ 
clude  the  following: 

a.  Isotherms  are  parallel  to  the  water  surface  both  laterally 
and  longitudinally. 

b.  The  water  in  each  discrete  layer  is  physically  homogeneous. 

c.  Internal  advection  and  heat  transfer  occur  only  in  the 
vertical  direction. 

d.  External  advection  (inflow  and  outflow)  occurs  as  a  uniform 
distribution  within  each  layer. 

e.  Internal  dispersion  (between  layers)  of  thermal  energy  is 
accomplished  by  a  diffusion  mechanism  that  combines  the 
effects  of  molecular  diffusion,  turbulent  diffusion,  and 
thermal  convection. 

4.  The  surface  heat  exchange,  internal  mixing,  and  advection  pro¬ 
cesses  are  simulated  separately,  and  their  effects  are  introduced  sequen¬ 
tially  at  each  time-step.  A  simplified  flow  chart  of  the  mathematical 


A3 


simulation  procedure  is  presented  in  Figure  A2.  Each  of  these  processes 
will  be  discussed. 


Figure  A2.  Simplified  WESTEX  flow  chart 
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Surface  Heat  Exchange 


5.  The  net  heat  exchange  at  the  surface  is  composed  of  seven  heat 
exchange  processes: 

a.  Shortwave  solar  radiation. 

b.  Reflected  shortwave  radiation. 

c .  Long-wave  atmospheric  radiation. 

d.  Reflected  long-wave  radiation. 

e.  Heat  transfer  due  to  conduction. 

f.  Back  radiation  from  the  water  surface. 

g.  Heat  loss  due  to  evaporation. 

6.  The  surface  heat  transfer  process  is  solved  in  the  WESTEX 
model  by  an  approach  developed  by  Edinger  and  Gpyer  (1965).  The  thermal 
equation  quantifying  the  net  surface  heat  exchange  (after  some  lineari¬ 
zation)  is: 


Hg  =  K(E  -  6g)  (A2) 

where 

2 

H  =  net  rate  of  surface  heat  transfer,  Btu/ft  /day 

a  2 
K  =  coefficient  of  surface  heat  exchange,  Btu/ft  /day/°F 

E  =  equilibrium  temperature,  °F 

0g  =  surface  temperature,  °F 

Equilibrium  temperature  is  defined  as  that  temperature  at  which  the  net 
rate  of  heat  exchange  between  the  water  surface  and  the  atmosphere  is 
zero.  The  coefficient  of  surface  heat  exchange  is  the  rate  of  heat 
transfer  at  the  air-water  interface.  The  computation  of  equilibrium 
temperature  and  heat  exchange  coefficient  is  based  solely  on  meteoroi '^i- 
cal  data  as  outlined  by  Edinger,  Duttweiler,  and  Geyer  (1968). 

7.  The  components  of  surface  heat  exchange,  with  the  exception  of 
shortwave  radiation,  are  immediately  absorbed  at  the  surface  or  in  the 
top  few  feet  of  the  pool.  Depending  upon  the  rolor  and  clarity  of  the 
water,  shortwave  radiation  penetrates  and  increases  the  temperature  at 
gre  ter  depths.  Based  on  laboratory  investigations,  Dake  and  Harleman 
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(1966)  have  proposed  an  exponential  decay  with  depth  for  describing  the 
heat  flux  due  to  shortwave  penetration.  This  approach  is  used  in  WESTEX. 

8.  The  surface  heat  exchanges  are  implemented  in  the  model  by 
the  placement  of  varying  percentages  of  the  incoming  shortwave  radiation 
in  each  layer  of  the  lake  and  by  the  placement  of  all  other  sources  of 
surface  heat  exchange  into  the  surface  layer.  The  shortwave  radiation 
is  distributed  exponentially,  so  that  most  is  absorbed  in  the  top  layers 
and  relatively  little  is  absorbed  in  the  lower  layers.  The  procedure 
can  be  expressed  mathematically  by  the  following  two  equations: 


Hg  =  K(E  -  8)  -  (1  -  m 


(A3) 


-A.z . 

H.  =  (1  -  p)tj»e  1  (A4) 

where 

2 

Hs  =  heat  >  msfer  rate  into  or  out  of  surface  lny«.r,  Btu/ft  /day 

P  =  perce  ge  of  incoming  solar  radiation  absorbed  in  surface 
layer 

*ji  =  total  incoming  shortwave  radiation  rate,  Btu/ft^/day 

2 

H.  =  rat>  of  heat  ausorpt  jn  in  layer  i  ,  Btu/ft  /day 

1  -1 
A.  =  light  extinction  coefficient,  ft 

z^  =  depth  of  layer  i  below  surface,  ft 


9.  Equations  A3  and  A4  are  applied  once  during  each  1-day  time- 
step.  The  net  heat  exchange  rate  into  each  layer  is  computed  and  con¬ 
verted  into  a  temperature  chan-e.  The  temperature  changes  are  used  to 
determine  an  updated  temperature  profile  for  the  lake. 


Inflow 


10.  The  inflow  process  is  simulated  numerically  in  three  basic 
steps.  The  point  of  neutral  buoyancy  of  the  inflow  is  found;  water  in 
the  lake  is  displaced  by  and  mixed  with  the  inflow;  and  a  new  water- 
surface  elevation  is  computed. 

11.  The  point  of  neutral  buoyancy  is  found  by  a  linear 
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interpolation  or  extrapolation  upon  the  density  profile  of  the  lake. 

The  inflow  volume  is  allocated  to  the  layer  of  neutral  buoyancy.  The 
contents  of  the  layer  of  neutral  buoyancy  are  then  fully  mixed  with  the 
inflow  quantity,  thereby  producing  a  volume-weighted  average  temperature 
for  this  layer.  If  the  inflow  volume  into  the  neutrally  buoyant  layer 
(layer  i)  causes  the  physical  capacity  of  that  layer  to  be  exceeded,  the 
excess  is  displaced  upward  at  the  mixed  temperature  of  the  inflow  layer. 
This  displacement  either  flows  into  the  next  higher  layer  (i+1)  or  forms 
a  new  surface  layer  (described  in  the  next  paragraph).  If  the  layer  of 
neutral  buoyancy  is  below  the  surface  layer,  the  excess  Is  fully  mixed 
with  layer  i+1  and  a  new  volume-weighted  temperature  for  that  layer  is 
produced.  This  process  continues  in  this  sequential  fashion  until  the 
introduction  of  the  excess  volume  from  one  layer  into  the  next  higher 
layer  does  not  exceed  the  physical  capacity  of-  the  upper  layer.  In  this 
manner  increments  of  inflow,  whose  magnitude  decreases  with  increasing 
distance  from  the  inflow  layer,  are  distributed  from  the  Inflow  layer  to 
the  surface. 

12.  If  the  inflow  current  is  found  to  be  an  overflow  (the  inflow 
density  is  less  than  that  of  the  surface  layer),  the  inflow  quantity  is 
mixed  with  the  volume  of  the  surface  layer.  If  the  Inflow  quantity  ex¬ 
ceeds  the  volume  of  the  surface  layer,  the  excess  forms  a  new  surface 
layer  at  the  mixed  temperature  of  the  inflow  layer.  The  addition  of  the 
inflow  quantity  in  any  manner  results  in  an  Increase  In  the  surface  ele¬ 
vation.  A  corresponding  decrease  in  the  surface  elevation  occurs  as  a 
result  of  the  outflow  simulation  process. 

Internal  Mixing 

13.  The  internal  mixing  process  is  represented  by  a  mixing  scheme 
mathematically  similar  to  a  simple  diffusion  analogy.  The  mixing  process 
is  used  to  account  for  the  effects  of  molecular  diffusion,  turbulent 
diffusion,  and  additional  internal  mixing  processes  which  are  not  ex¬ 
plicitly  addressed.  Internal  mixing  transfers  heat  and  other  water  qual¬ 
ity  constituents  between  adjacent  layers.  The  magnitude  of  the  transfer 
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between  two  layers  in  the  WESTEX  model  is  expressed  as  a  percentage  of 
the  total  transfer  required  to  mix  the  two  layers.  This  percentage  is  a 
mixing  coefficient  which  is  defined  for  every  layer.  Data  input  includes 
values  of  the  mixing  coefficient  at  the  maximum  pool  elevation  and  at 
the  bottom  of  the  lake.  The  values  are  predicated  upon  field  data.  An 
exponential  fit  between  the  two  extreme  values  is  used  to  determine  the 
appropriate  coefficient  for  each  layer.  The  same  coefficients  are  used 
for  mixing  other  water  quality  constituents  as  for  mixing  heat. 

Outflow 


14.  The  outflow  component  of  the  model  incorporates  the  selective 
withdrawal  techniques  for  orifice  flow  developed  at  WES  by  Bohan  and 
Grace  (1973).  Transcendental  equations  defining  the  location  of  the 
zero-velocity  limits  are  solved  iteratively.  The  zero-velocity  limits 
are  functionally  dependent  on  the  release  flow  rate  and  the  in-lake  den¬ 
sity  structure.  With  knowledge  of  the  withdrawal  limits,  the  outflow 
velocity  profile  can  be  determined,  and  the  change  in  the  internal  heat 
budget  can  be  quantified. 

15.  The  relationship  used  to  define  withdrawal  limits  for  flow 
through  orifices  is  the  following: 


Q  = 


(A5) 


where 

Q  =  total  flow  through  the  port,  cfs 

Z  =  vertical  distance  from  the  orifice  center-line  elevation  to 
w 

withdrawal  limit,  ft 

Ap  =  density  difference  between  elevation  of  orifice  center  line 
and  withdrawal  limit,  g/m£ 

p  =  fluid  density  at  orifice  center  line,  g/m£ 

0  2 
g  =  gravitational  acceleration,  32.18  ft/sec 

It  is  not  possible  to  solve  Equation  A5  in  closed  form  for  Z^  ,  but  it 

can  be  solved  numerically. 
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16.  After  determination  of  upper  and  lower  withdrawal  limits,  the 
location  of  maximum  velocity  is  determined  from  the  following  equation: 


H 


w 


(A6) 


where 

Yj  =  distance  from  lower  limit  to  location  of  maximum  velocity,  ft 

=  thickness  of  withdrawal  zone;  distance  from  lower  limit  to 
upper  limit,  ft 

=  distance  from  lower  limit  to  orifice  center  line,  ft 

The  velocity  profile  above  or  below  the  elevation  of  maximum  velocity 
can  be  calculated  for  boundary  interference  by: 


V 

V 

max 


1 


(A7) 


and  for  the  case  when  the  withdrawal  zone  is  not  limited  by  a  boundary: 


V 

V 

max 


(A8) 


where 

V 

xV 

max 

y 


Y 


Ap 


Ap 


max 


ratio  of  local  velocity  to  maximum  velocity 

distance  from  elevation  of  maximum  velocity  to  local 
elevation,  ft 

distance  from  elevation  of  maximum  velocity  to  the  with¬ 
drawal  limit  of  interest,  ft 

density  difference  between  elevation  of  maximum  velocity 
and  local  elevation,  g/m£ 

density  difference  between  elevation  of  maximum  velocity 
and  the  withdrawal  limit  of  interest,  g/m2 


If  multiple-level  ports  are  open,  then  a  flow-weighted  relative  velocity 


profile  is  computed  independently  for  each  port,  and  the  velocity  pro 


files  are  superimposed  on  the  basis  of  a  controlled  shift  of  the 
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withdrawal  limits  in  the  zone  of  overlap  to  achieve  a  total  relative 
velocity  profile. 

17.  Lake  regulation  algorithms  have  been  developed*  to  realisti¬ 
cally  simulate  the  field  operation  of  a  selective  withdrawal  system. 

The  selective  withdrawal  system  is  assumed  to  be  configured  with  an 
arbitrary  number  of  selective  withdrawal  intakes  located  in  each  of  two 
wet  wells  with  a  separate  floodgate.  Maximum  flows  and  minimum  flows 
from  each  intake  and  from  the  floodgate  must  be  specified.  Also,  the 
maximum  flow  for  the  selective  withdrawal  system  is  specified.  The  algo¬ 
rithms  attempt  to  numerically  withdraw  water  at  or  near  the  objective 
temperature.  Withdrawal  will  be  from  either  one  intake  level  or  two  ad¬ 
jacent  intake  levels,  and/or  the  flood-control  intake  depending  upon  the 
objective  temperature,  the  temperature  profile,  the  intake  capacities, 
and  the  amount  of  flow  to  be  released. 

Density  Stability 

18.  Cooling  of  the  lake  surface  causes  a  density  instability  and 
results  in  convective  mixing  within  the  water  column.  Stability  is 
checked  by  searching  adjacent  layers  from  bottom  to  top  and  comparing 
densities.  If  a  density  instability  is  identified,  the  two  unstable 
layers  are  mixed,  and  the  mixed  density  is  compared  with  the  density  of 
the  layer  above  the  mixed  region.  If  an  instability  still  remains,  the 
layer  above  the  mixed  region  is  included  in  the  mixed  region,  and  the 
process  continues  until  stability  is  achieved  or  the  surface  is  reached. 
By  mixing  layers  above  an  instability  it  is  possible  to  create  an  insta¬ 
bility  below  the  mixed  region.  If  such  an  instability  is  detected,  then 
mixing  proceeds  downward  until  stability  is  achieved  or  the  bottom  is 
reached. 


*  See  footnote  reference,  page  9,  main  text. 
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APPENDIX  B:  FORMULATION  OF  OPTIMIZATION  METHODOLOGY 


Formulation  of  Optimization  Procedure 

1.  Initial  simulation  of  Cowanesque  Lake  with  increased  normal 
pool  elevations  indicated  a  need  for  the  location  of  additional  selective 
withdrawal  intakes  to  meet  downstream  objective  temperature  objectives. 
The  location  of  these  additions  to  the  present  selective  withdrawal  sys¬ 
tem  (which  houses  two  ports  at  el  1037.0  and  two  at  el  1014.5)  was  ac¬ 
complished  by  the  following  sequential  process: 

£.  Estimate  the  location(s)  of  the  additional  port(s)  needed 
for  minimization  of  downstream  temperature  violations. 

1).  Combine  the  estimated  port  location(s)  with  the  existing 
system.  This  intake  configuration  was  used  as  an  input 
into  the  WESTEX  thermal  model  to  simulate  the  in-lake  and 
release  characteristics  resulting  from  operation  of  this 
intake  scheme. 

£.  Compute  an  objective  function  (described  in  detail  in 

paragraph  4)  based  on  deviation  of  objective  temperatures 
and  release  temperatures  of  step  t>. 

d_.  Predict  a  new  location  estimate  for  the  additional  port(s) 
based  on  steps  _a,  _b,  and  c_. 

Steps  J),  £,  and  id  are  repeated  until  a  minimum  value  for  the  objective 
function  of  step  £  is  obtained.  At  the  minimum  objective  function  the 
location  of  the  additional  selective  withdrawal  intakes  is  the  "optimum" 
location  for  maintenance  of  downstream  temperature  objectives. 

2.  Although  this  procedure  could  be  accomplished  manually,  tech¬ 
niques  of  mathematical  optimization  systematically  examine  the  possible 
estimates  for  port  location  and,  in  general,  converge  to  the  best  de¬ 
cision  with  the  least  expense  of  computer  resources.  In  this  study  an 
optimization  routine,  requiring  the  output  of  steps  £,  _b,  and  £  as  input, 
was  used  to  locate  an  estimate  to  the  "optimum"  selective  withdrawal  con¬ 
figuration  for  each  operating  condition  (step  d) .  The  optimization  pro¬ 
cedure  continues  to  iterate  until  the  difference  between  the  elevations 
of  the  updated  intake  configuration  and  the  corresponding  elevations  of 
the  intake  configuration  at  the  previous  iteration  is  less  than  a 
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predefined  tolerance  (0.25  ft  for  this  study).  The  objective  function 
is  then  minimized  for  the  given  operating  condition  and  the  elevations 
of  the  intake  configuration  from  the  previous  iteration  are  defined  to 
be  the  "optimum.'' 


Optimizer  Description 

3.  The  optimization  routine  used  in  this  study  finds  the  minimum 
of  a  general  function  f(x)  of  a  single  variable  without  using  deriva¬ 
tive  values.  The  routine,  NW034  from  the  Boeing  Computer  Services  Soft¬ 
ware  library,  minimized  an  objective  temperature  function  whose  indepen¬ 
dent  variable  was  the  elevation  of  a  single  level  of  additional  selective 
withdrawal  intakes.  Following  initialization,  the  procedure  requires 
knowledge  of  three  objective  function  evaluations  and  the  points  (port 
locations)  at  which  they  are  evaluated.  Using  parabolic  interpolation, 
NW034  finds  a  local  minimum  of  a  parabola  fitted  to  these  three  sets  of 
points.  The  data  pair  (elevation,  function  value)  having  the  maximum 
function  value  is  then  discarded.  A  new  parabola  is  developed  from  the 
two  remaining  pairs  plus  a  new  data  pair  consisting  of  the  local  minimum 
of  the  previous  parabola  and  the  objective  function  evaluated  at  this 
minimum.  This  procedure  continues  until  the  elevation  difference  be¬ 
tween  two  successive  local  minimums  is  less  than  some  tolerance  (con¬ 
vergence)  or  a  maximum  number  of  iterations  has  been  reached.  The 
"optimum"  additional  port  level  location  is  the  elevation  value  at  the 
local  minimum  found  for  tne  iteration  previous  to  convergence. 

Development  of  Objective  Function 

4.  The  principal  purpose  of  this  study  was  to  determine  location 
of  additional  selective  withdrawal  intakes  to  meet  a  specified  down¬ 
stream  warrawater  temperature  objective.  An  objective  function  is  a 
scalar  index  that  provides  an  indication  of  how  good  one  possible  deci¬ 
sion  (i.e.,  port  location  estimate)  is.  Minimization  of  this  objective 
function  produces  the  optimal  locations  of  these  additional  ports. 
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5.  The  objective  function  used  for  Cowanesque  Lake  is  the  weighted 
sura  of  several  terras.  The  principal  terms  are  weighted  portions  of  the 
sum  of  squares  of  deviations  between  computed  release  and  objective  tem¬ 
perature.  The  primary  terra  represents  the  contribution  of  the  squared 
temperature  deviation  on  days  when  the  absolute  value  of  the  deviation 
(release  temperature  from  objective  temperature)  is  greater  than  2.78°C. 
This  portion  has  physical  significance  in  that  the  State  of  Pennsylvania 
requires  releases  from  Cowanesque  Lake  to  be  plus  or  minus  2.78°C  from 
the  natural  stream  temperature  (which  is  the  study  objective).  Thus,  on 
days  when  the  release  temperature  is  within  this  band,  no  contribution 

to  the  sum  of  squares  is  recorded.  This  terra  was  found  to  minimize  devi¬ 
ation  over  the  entire  simulation  period  much  better  than  the  total 
squared  deviation  between  the  release  and  objective  temperatures  summed 
over  each  day  of  the  period. 

6.  A  secondary  term  represents  the  contribution  of  the  total  sum 
of  squared  deviations  between  release  and  objective  temperature.  Several 
intake  configurations  were  found  that  minimized  deviation  of  the  release 
temperature  from  the  ±2.78°C  band.  This  secondary  term  was  used  to  give 
mathematical  uniqueness  to  each  of  these  configurations  by  adding  a  con¬ 
tribution  to  the  objective  function  for  the  deviation  of  the  release  tem¬ 
perature  from  the  objective  itsPif,  and  not  merely  the  band,  for  every 
day. 

7.  Two  further  terms  were  added  to  the  objective  function  to  re¬ 
strict  the  range  of  possible  decisions.  Each  of  the  terms  penalized  the 
objective  function  for  certain  port  placements.  The  first  term  penalized 
the  objective  function  for  placement  of  a  port  whose  top  elevation  was 
above  the  water  surface.  The  second  term  added  a  penalty  for  port  place¬ 
ment  below  el  1037.0  because  tne  present  system  of  selective  withdrawal 
ports  at  Cowanesque  is  adequate  to  handle  withdrawal  at  this  elevation 
and  below. 
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